Broadband Mueller ellipsometer as an all-in-one tool for spectral and temporal analysis of mutarotation kinetics

Spectroscopic Mueller matrix ellipsometry is becoming increasingly routine across physical branches of science, even outside optics. The highly sensitive tracking of the polarization-related physical properties offers a reliable and non-destructive analysis of virtually any sample at hand. If coupled with a physical model, it is impeccable in performance and irreplaceable in versatility. Nonetheless, this method is rarely adopted interdisciplinarily, and when it is, it often plays a supporting role, which does not take benefit of its full potential. To bridge this gap, we present Mueller matrix ellipsometry in the context of chiroptical spectroscopy. In this work, we utilize a commercial broadband Mueller ellipsometer to analyze the optical activity of a saccharides solution. We verify the correctness of the method in the first place by studying the well-known rotatory power of glucose, fructose, and sucrose. By employing a physically meaningful dispersion model, we obtain 2π-unwrapped absolute specific rotations. Besides that, we demonstrate the capability of tracing the glucose mutarotation kinetics from just one set of measurements. Coupling the Mueller matrix ellipsometry with the proposed dispersion model ultimately leads to the precisely determined mutarotation rate constants and spectrally and temporally resolved gyration tensor of individual glucose anomers. In this view, Mueller matrix ellipsometry may stand as an offbeat yet equal technique to those considered classical chiroptical spectroscopy techniques, which may help open new opportunities for broader polarimetric applications in biomedicine and chemistry.


Introduction
In recent years, cutting-edge applications of Mueller matrix spectroscopic ellipsometry (MMSE) have been bursting far beyond its homeland of photonics and material sciences, bringing us probes into cancer diagnostics, 1-3 ultrafast electronic processes, 4,5 in vivo imaging, 6,7 solar and hydrogen energy harvesting, [8][9][10] toxin detection, 11 and much more. The unique exibility of MMSE stems from the contactless ability of light polarization to track even the subtlest material attributes. It is virtually the signature of the interplay between circular polarizations and chiral objects, providing us with experimental observables and elucidating their spatial conformation at the molecular level. Specically, different travel speeds of le-and right-circular polarizations in a given medium reect its optical activity (OA). In the light of the MMSE, we directly observe a Mueller matrixa characteristic quantity embracing OA dispersion (optical rotation) and its Kramers-Kronig related absorptive counterpart (circular dichroism).
The leading chiroptical techniques are mainly vibrational circular dichroism 12,13 (VCD), Raman optical activity [14][15][16] (ROA), and electronic circular dichroism 17 (ECD), which, in conjunction with ab initio calculations, 18 reveal the absolute conguration of a molecule. Within those, it seems quite ambitious to justly pick MMSE and challenge the well-established methods. However, contemporary ellipsometers have come a long way from the traditional single-wavelength polarimetry (saccharimetry), particularly in greater spectral range, high accuracy, capability of real-time tracking of the optical response, efficient data processing, and much more. Bottom-line, putting the polarimetric techniques into the context of state-of-the-art MMSE would open up new opportunities in chiroptical spectroscopy.
Covering the frequency range typically from ultraviolet to the infrared spectrum, the Mueller matrix reects both electronic and vibrational transitions responsible for the chiroptical response of a chiral molecule. As a rule of thumb, the ellipsometric community has adopted the glucose molecule as a universal representation of a chiral object within the landscape of chemistry. Its clinical implications have attracted attention exploiting MMSE versatility by sensing glucose 19 in  [20][21][22] body-uid environments, 23 or even phantom tissues 24 and actual human tissues. 25 Moreover, fast and precise response makes MMSE widely popular as an in situ control 26,27 capable of real-time tracking of dynamic processes. This has been done very sparsely in chemical kinetics, 28 and only once regarding glucose mutarotation, 29 despite being fundamental to saccharides. Considering that the used method is monochromatic and does not account for glucose gyration spectra, it leaves much room for improvement, where MMSE appears to be a hot candidate.
In this article, we present commercially available MMSE as an up-to-date tool for non-destructive sensing of saccharide chiroptical response capable of tracking the mutarotation kinetics. First, we will show the power of the spectral Mueller matrix in tandem with a physically consistent model. The correct t to the experimental data enables retrieval of the isotropic gyration tensor components spectra, gaining the sensitivity to 2p-unwrapped absolute optical rotation. A simple combination of the spectral approach with the rate equation promotes the proposed method to the next level. We will show that the MMSE tracing of just a single forward mutarotation reaction is sufficient enough to reveal (i) absolute optical rotations at all measured wavelengths at any time instant, (ii) spectrally and temporally resolved effective gyration response of the measured solution, (iii) mutarotation rate constants and corresponding concentrations evolution, and, remarkably, (iv) gyration spectra of individual aand b-anomers of D-glucose. These demonstrate the paramount strength of the proposed method as a chiroptical technique for quick, easy-to-use, and reproducible analysis of various static and dynamic processes in physical chemistry, potentially providing an experimental counterpart to specic quantum-mechanical quantities employed in ab initio modeling of OA. 30 Note, right before we focus on the main results (Section 3), the next section introduces the necessary technicalities of the method.
Before measuring the spectral responses of D-glucose and Dfructose, their solutions were rst le at rest for 24 hours to reach the state of their anomeric equilibrium at a given temperature. Solutions of D-sucrose, which is a non-reducing and therefore mutarotation inactive saccharide, were measured just aer 4 hours, only to ensure perfect dissolution. Mutarotation ellipsometric experiments on a-D-glucose were performed right away with freshly prepared solution.

Mueller matrix ellipsometry
The method principle reects its very namethe measurement of the change in the arbitrary (generally elliptical) polarization state of light aer interaction with a sample. The polarization state change records a 4 × 4 Mueller matrix M, which reects the complete polarimetric response of a sample, including the effect of optical activity in chiral liquids.
The corresponding reduced (normalized to m 11 ) Mueller matrix M M ¼ (1) accounts for the inuence of OA in the form of its dispersive part G(l) and absorptive part G ′ (l). The OA spectral (wavelength, l) dependence is connected with the different propagation speeds of le (L) and right (R) circularly polarized light, or eigenmodes (those states of a polarized light that do not change their polarization state when passing through a given environment), in a medium of length d. It is caused by different dispersion of L and R eigenmodes refractive indices n L,R . The delay between the eigenmodes, called circular phase retardation, effectively manifests itself as a rotation of the plane of linearly polarized light (see Fig. 1a) over the angle which corresponds to the traditional denition of OA. Refer to Appendix A.1 for details. Of course, this denition is consistent with the theory on the microscopic optical activity approach accounting for chiral molecules with random orientation dissolved in an achiral solution, refer to Appendix A.2 for details. All experimental work was done on a standard Woollam RC2-DI (dual-rotating compensator) Mueller matrix ellipsometer in transmission conguration with the spectral range of 193 nm-1700 nm (6.42 eV-0.73 eV). Its instrumentation is based on a polarization state generator (PSGprobing arm: source, polarizer, and rotating compensator) and polarization state analyzer (PSAdetection arm: rotating compensator, analyzer, detector part).

Temperature cell design and control
For temperature-controlled measurements, we have designed a custom Peltier device-based 3D-printed cuvette holder. Fig. 1b shows the cell schematics. The keystone element of the temperature control is a pair of Peltier modules (HP-127100HTS, 90 W, temperature range from −20°C to +200°C) with each module attached to one side of the optical glass cuvette with the optical path length 50.04 mm working as a heat pump. The Peltier modules are powered by a current source (BOP20-10M Kepco Bipolar Power Supply) in a parallel circuit. The actual temperature of the measured solution is sensed in real-time by a PT100 temperature sensor connected to a PT-104 data logger.
We connected the outer sides of the Peltier modules to a pair of Corsair H45 computer water cooling units to ensure their stable temperature. The second critical aspect of our setup is a good thermal conduction between each contact of the Cuvette-Peltier-Cooler unit. To do so, we encased all these parts in a 3D-printed exoskeleton allowing us to bolt them together and apply a thermal paste. The system is equipped with a small magnet attached to an electric motor under the cuvette, which can be used as a magnetic stirrer to ensure a homogeneous temperature distribution in the solution.

Mutarotation experiment methodology and rate equation
Mutarotation is a change in the concentrations of the a-anomer and b-anomer of reducing saccharides upon dissolution in water. The presence of the hemiacetal carbon causes, through ring-chain tautomerism, a gradual inter-conversion between the aand b-anomers until equilibrium is established. It results in a change in the optical rotation of a given solution. Fig. 3a shows the schematics of the glucose mutarotation. It is not oen pointed out that a tiny fraction of the six-membered rings (pyranoses) convert to ve-membered rings (furanoses) during mutarotation, regardless of the starting substance from which the solution is prepared. The system at equilibrium remains a dynamic environment in which inter-conversions are constantly occurring over and over, yet the overall concentration ratio among all forms is conserved. The rst effect mentioned is negligible as far as the concentrations of the ve-membered rings are concerned, and the second doesn't affect the experiment. Therefore, we do not include either effect in the description of mutarotation.
The mutarotation of D-glucose is described by the reverse reaction where s 0 and s * are forward and reverse mutarotation rate constants, respectively. The reaction kinetics evaluates the change in aand b-anomer concentrations C a,b and is given by a rst-order rate equation with initial conditions and solution where C a,eq is the a-anomer concentration in equilibrium. The observed total rate constant characterizes the time evolution of D-glucose aand b-forms instantaneous concentrations C a,b via forward and reverse rate constants s 0 , s * , respectively, and equilibrium constant K,  (2) with adjustable input current (3), which provide temperature stability inside the cuvette with the measured solution (4). The temperature is sensed by a PT100 sensor (5) and recorded in realtime (6).  We set the initial time t = 0 of the a-D-glucose mutarotation process to the instant at which crystals were no longer visible in the aqueous solution. It takes approximately 30 s for concentrations of C a = C = 0.25 M. The rst ellipsometric measurement was taken aer 300 s. We monitored the temperature T inside the cuvette throughout the process with the accuracy of 0.5°C for the following 12 hours, at which point the experiment was terminated.  (2) does not constitute an optimal model for subsequent minimization against the experimental data, since the difference n L − n R incorporates the effective cumulative response of both medium refractive index dispersion n and chirality dispersion.

Spectral model and method validity
Therefore, we propose to apply the dispersion model that combines an optics-based approach based on oscillator models, 34,35 accounts for material symmetry dening the shape of material (permittivity, gyration) tensors, 36,37 and employs rigorous wave-equation-based calculations of the polarization eigenmodes, 38,39 The model separates the gyration tensor element g 11 from the medium refraction index n, both described by phenomenological constants A, B, A i , B i , of which the rst two had set as the tting parameters. The constants A i , B i describing n of water were taken from ref. 40. Note that the refractive index of water differs from that of a mixture of water and saccharide, however, the difference is negligible for dilute solutions 41 -5% change in saccharide to water mass ratio causes less than 1% change in the refractive index dispersion (the mass ratio of 0.25 M saccharide solution is approximately 4.5%). The spectral limitation of the measurement comes from the transparency of the optical glass cuvette (>300 nm) and the water absorption (>1000 nm). Since monosaccharides exhibit electronic circular dichroism in the spectral region below 190 nm, this is a phenomenon outside the spectral range of the ellipsometer employed. Thus, G ′ (l) = 0. The striking advantage of having spectral measurements is the sensitivity to the absolute (2kp-wise, k˛ℤ ) optical rotation of studied media. It gives the exibility to measure very high concentrations of the solution (several mol l −1 ), compounds with enormous optical activity (e.g., helicenes), or both. Fig. 2a shows a practical example. In the spectra of M of the 3 M solution of D-fructose, we have marked points P1 and P2 corresponding to measurements at randomly picked wavelengths of 382 nm, and 627 nm. The rst advantageous feature of the spectral Mueller matrix M is the complementarity between the elements m 22 , m 33 and m 23 , m 32 resulting from the evenness and oddness of the sine and cosine functions. Thus, based on the combination of the signs of the analyzed elements of M, it is possible to unambiguously determine G(P1, P2) corresponding to the spectral points P1 and P2 in Fig. 2a, yet only within the 2p range. The second aspect is the dispersion model. Its parameters are specic for every material, varying in the amplitude, curvature, and offset of the observed dispersion, and their combination is unique to each point in the spectrum. These In addition, the dispersion model enables us to distinguish the data from noise or systematic measurement errors even for particularly weak signals. Considering that the Woollam RC2 theoretical threshold value for distinguishing signal from noise in the Mueller matrix is 0.001, the lowest detectable G = 0.5 mdeg. This value corresponds to the lowest measurable concentrations of D-glucose (19.0 mg ml −1 ), D-fructose (10.8 mg ml −1 ) and D-sucrose (15.3 mg ml −1 ) at 589 nm. Fig. 2b shows D-glucose, D-fructose and D-sucrose optical rotatory power spectra [G] T l dened as We have particularly indicated [G] 20 589 (measured at the temperature T = 20°C and l = 589 nm), as these are generally wellknown values of specic rotatory power obtained by monochromatic setups working with a sodium D-line lamp. A comparison of the measured values with those reported in the literature is provided in Table 1 and shows quite a good agreement between them. The deviations could be mainly due to minor errors caused by weighing and solution preparation. Fig. 2c shows the tted g 11 dispersions of the studied saccharides.
Up until now, we have been considering the effective cumulative optical response and describing the studied materials accordingly. In the above-studied cases, the physical meaning of corresponding dispersions g 11 (see Fig. 2c) is not directly the response of pure substances (pure anomers) but rather an effective value for a given solution in equilibrium, which may consist of an arbitrarily high number of chiral chemical individuals. Given this logic, each individual in the solution manifests its own g 11 . This approach is crucial for the study of mutarotation, as it gives us access to all anomers involved.

Application of the ellipsometry to mutarotation kinetics
We assume that the observed gyration response of the solution is proportional to the concentration-weighted gyration response of the individual anomers in the solution: Gðl; tÞ ¼ p lnðlÞ where g 11,a and g 11,b are gyration tensor dispersions of pure aanomer and b-anomer of D-glucose. Such a combination of the spectral model (9) and the rate equation (6) allows for a physical description of G(l,t) based on (2). 3.2.1 Spectral dependence. Fig. 3b shows model (eqn (11), red curves) optimized to the experimental data (blue curves) of the a-D-glucose mutarotation process. Each curve stands for the [G] 20 l time evolution of 0.25 M a-D-glucose solution during ongoing mutarotation. For better clarity, we show results for selected wavelengths onlythe actual measured and modeled spectral range is 300-1000 nm with the step of 1 nm. As expected, the observed spectra follow the exponential tendency given by (6) and converge to the value of [G] 20 l (t / N) corresponding to the equilibrium state between the anomers in the solution. The inset of Fig. 3b shows the [G] 20 589 (t) curve converging to the value [G] 20 589 (t = 720 min) z 52 deg ml g −1 dm −1 , which is consistent with those in Table 1. Now, with the modeled spectra available, we dig more into their physical signicance. We present the individual gyration tensor contributions g 11,a (l) and g 11,b (l), see Fig. 3c. In the following, we will address the validity of our ndings: the ellipsometric measurements start with a freshly prepared solution of a-D-glucose and the measured data initially contains mainly information about the dispersion of g 11,a , which we are thus able to model accurately. As the mutarotation reaction progresses, the effect of g 11,b dispersion gets more dominant in the optical response of the system ever closer to its equilibrium state. The consequence is that the g 11,b is not accessible in a portion of measurements causing lower numerical sensitivity to this contribution, which accounts for the correlation between the parameters describing g 11,b and C a,eq . In order to correctly determine the dispersion of g 11,b at all, we have to assume the known value of C a,eq , in this case C eq,a : C eq,b = 36% : 64% at T = 20°C. Elsewise, the dispersions of g 11,b and C a,eq correlate in the model. This step bears the risk of possible inaccuracies and should be justied: applying the discussed ratio, the gyration response of the equilibrium solution is equal to the concentration-weighted sum of the gyration spectra of the isolated anomers (blue and red curves in Fig. 3c) resulting in the yellow curve. In other words, this curve must reproduce (and it actually does with 97.5-99% accuracy) the gyration response of previously measured D-glucose (black dotted curve) as in Fig. 2c 20 589,b = 18.7 deg ml g −1 dm −1 . Clearly, there are some reservations in the calculated absolute values, but they may stand for a good approximation placing the results into the frame of what to expect. In any case, the extracted spectral dependence g 11,a,b (l) may helpfully stand as the experimental counterpart to the quantities proportional to electric dipoleelectric dipole and electric dipole-electric quadrupole optical activity tensors used in theoretical simulations. 30 3.2.2 Temporal dependence. Having established the spectral response throughout the mutarotation, we now focus on modeling the kinetics. Fig. 3d provides a graphical representation of the C a , C b time evolution, as the mutarotation conversion between the D-glucose anomers proceeds. The rate constants s, s 0 , s * characterize the speed of this conversion. The s is a directly tted parameter of model (11), in which it characterizes the exponential decay (6)  As it is not possible to directly compare the latter two methods with our approach, we do not include them in our further discussion.
The deviations from the tabulated values could come from several causes. Trivially, inevitable errors are due to possible inaccuracy in the concentration of the initial solution, and temperature dependence of G and [G] 20 l . Moreover, the model does not account for the effect of the impurities in the initial powders, 96% a-D-glucose in particular. We adopt a standard approachthe composition of impurities is expected to be a racemic mixture of D-glucose, which certainly has a minor contribution to the data at the beginning of the experiment. In the view of the method demonstration and to keep the model reasonably simple, we omit this. The crucial part rests in the accuracy of s, which depends on the accuracy of G (G prole forms the mutarotation exponential). Various authors report G precision of 8.3 × 10 −5 deg cm −1 , 29 and 0.01 deg per 40 cm. 42 The heterodyne polarimeter used in ref. 29 is more accurate than our ellipsometer with an accuracy of 5 × 10 −4 deg per 5 cm. Nevertheless, at this scale, it doesn't cause such differences among reported rate constants. Even the relatively strong deviations in G(t) (forming the shape of the mutarotation exponential point-by-point {t 1 , t 2 , .} and thus introducing the possible errors in a pointwise manner) cause a relatively small inaccuracy in s characterizing the decay of the whole exponential formed as the G(t i ). Nonetheless, the coupling between spectral approach and rate equation resulting in the model (11) is what gains the upper hand over other, usually single-wavelength (typically 589 nm and 632.8 nm), polarimetric approaches. We t the s simultaneously over the hundreds of exponentials, each for a different wavelength (see Fig. 3b). In other words, we calculate s from a much bigger dataset than is usually considered in other studies such as those in Table 2. Moreover, G(l) follows a classical optical dispersion tendency, which results in stronger effects over the ultraviolet spectral region and weaker effects over the infrared. It benets the signal accessibility over shorter wavelengths, particularly for rarely considered l < 589 nm, providing experimental data less burdened by noise and systematic measurement errors. Monochromatic methods do not allow that at all. Thus, the robustness of our method signicantly enhances the precision and sensitivity with which we calculate the rate constant s.

Discussion and conclusion
In this work, we have presented an oeat implementation of custom extended spectral Mueller matrix ellipsometry to describe the chiroptical response of simple saccharides in conjunction with evaluating corresponding mutarotation kinetics. We chose glucose, fructose, and sucrose as the method demonstration structures. We showed that the measurements resulted in a characteristic form of the Mueller matrix whose optical dispersion we have described by a physically consistent model coupled to a rate equation. We show that the spectraltemporal coupling in the model enhances the reliability, precision, and interpretation of the revealed spectra, thus enabling accurate evaluation of both the essential polarimetric quantities and the mutarotation rate. Remarkably, we extracted the gyration response of both D-glucose anomers from a single forward mutarotation reaction.
Above all, the presented method represents a universal approach applicable beyond the compounds in this study. An initial suggestion would be to adopt this method for compounds with a similar gyration response, G ′ (l) = 0, such as more complex saccharides. In this case, our method is ready-to-use straightaway. In principle, it would be possible to determine the individual optical response of a solution of more than one saccharide at their unknown concentrations. In such a case, however, prior knowledge of the gyration response of the individual compounds of the mixture would be necessary due to the existence of multiple combinations of tted parameter values that would (numerically) equally describe the measured datait is clear, that the method presented cannot replace qualitative analyses such as chromatography, but may stand as their complementary tool. In particular, it is worth extending the measurements to substances with G ′ (l) s 0. In the spectral region given by the ellipsometer, these would typically be substances with an active chromophore and ideally exhibiting exciton coupling, enhancing the strength of the observed effect in the Mueller matrix elements m 14 and m 41 . The spectra (12) would then be Kramers-Kronig consistent. In contrast, measurements involving substances with a hierarchical arrangement (proteins) are technically measurable, but the presented model would be inadequate to describe them. It is due to the extra contribution to overall optical activity from its secondary (tertiary) structure. Still, the effective gyration response would remain accessible, which may keep the approach more than sufficient for lots of applications. Of course, the shape of the Mueller matrix (15) would be preserved in this case as well.
We believe that despite some technical drawbacks arising from the need to navigate within the Mueller calculus framework, our approach could stand as a quick, accurate, and repeatable tool pushing the bounds of traditional polarimetry well beyond its limitations. The Mueller matrix allows simultaneous access to G(l) and G ′ (l) separately from each other and even from other polarization effects. This paramount strength of the Mueller matrix approach has the inherent potential to open up a novel avenue for data interpretation and validation in conjunction with well-established techniques encompassing ROA, VCD, ECD, and ab initio calculations. For example, accurate evaluation of chemical kinetics may pave the way toward comprehensive polarimetric models for tissue characterization.

Appendices
A An explanatory insight into the physical aspects of the method We are aware of the possible discouragement of the reader from applying this method in practice due to the seemingly complexlooking theory surrounding Mueller matrices. Therefore, in this section, we attempt to present the essential polarimetric principles of the present approach in an accessible and pictorial way to the reader who may not be an expert in polarimetry and ellipsometry.
A.1 Mueller matrix of optically active liquids. Le and right circularly polarized eigenmodes (polarization states that do not change upon interaction with given environment) propagate in a chiral liquid medium of length d in the form of a plane monochromatic (single-wavelength, l) wave, characterized by wavenumber k 0 = 2p/l, with phase factor exp(−ik 0 ñ L,R d). Here, refractive indices n L,R and absorption coefficients n 0 L;R are the respective real and imaginary parts of the complex propagation constantñ L;R ¼ n L;R À in 0 L;R . The difference between n L,R effectively causes an azimuthal rotation of linearly polarized light plane about angle G, while G ′ is an angular measure of the difference between n 0 L;R (or the ellipticity angle) related to the eigenmodes absorption in the medium, The polarimetric response of the studied medium is given by transmission Jones matrix J. It might be hard to navigate within the Jones and Mueller matrix interpretation, if the medium possesses a complex polarimetric response. In specic casestransmission conguration, the observed effects are smallwe can think of the complex polarimetric response as a cascade of individual fundamental effects, as if the medium consisted of multiple pieces in a series, each with its own polarimetric effect. The corresponding Jones matrices multiply. For our case, this is represented in Fig. 4a we introduce separately the effects of absorption (diattenuation, D) and dispersion (retardation, R) of circularly polarized light. The matrix J LR (transmission Jones matrix of optically active medium expressed on the basis of circularly polarized eigenmodes) then takes the form of a product Eqn (13a) is valid, because the matrices J R , J D commute. Now, by transforming the J LR matrix into its Cartesian form J, applying some basic algebra and trigonometric identities (for detailed derivation, see ref. 46), one can obtain the following expression, Following the well-known procedure of Jones-to-Mueller : (15) In non-absorbing spectral regions, (G ′ = 0), the matrix (14) transforms into a matrix representing standard Cartesian rotator J = J R , and the corresponding Mueller matrix takes the form If the circular dichroism is small, it allows us to perform a Taylor expansion around the zero point, implying : A.2 Chiral molecules in isotropic solution. Here we address a potential concern that may occur to some as a result of a fairly intuitive notion: since the orientation of all polarizable chiral molecules in solution is random, shouldn't this lead to all temporal polarizations of chiral molecules averaged to zero in a given direction? If this were true, then ellipsometric measurements would be entirely inadequate for our application and would give misleading interpretation of G. To put this false doubt to rest, we offer the following explanation.
It relates the microscopic molecular quantities (in a 123 coordination system) with their macroscopic observables (in the xyz reference laboratory frame). To do so, we need to account for the various multipole moments within the ensemble of N molecules induced by the electromagnetic wave (characterized by the electric intensity E and magnetic ux density B) propagating in a solution. In the case of optically active liquid, the electric dipole m and magnetic dipole m approximation is sufficient, and the respective moments take the form 48-50 Fig. 4 (a) The total optical response of the chiral solution can be thought of as a commuting product of the empty cuvette (linear birefringence caused by the residual stress in the glass, LB), the optically active solution (the circular birefringence causing phase retardation G, CB), and the dichroic solution (the differential absorption of left and right circularly polarized light, G ′ , CD). Please take note that this is a great way for experimentalists to understand the experimental data, however, it must be used very carefully when it comes to the rigorous interpretation of the data. The observed effects must be small and measured in transmission configuration. Each effect is recorded separately within a different Mueller matrix element. (b) The graphical illustration of the Mueller matrix shows the interpretation of the relevant matrix elements (main diagonal disregarded). While the minor diagonal elements always show a chiroptical response, the green LB terms would drop to zero when using stress-free cuvettes.
where the molecular polarizabilityâ and optical activity tensorb are in the general form Since the orientation of each molecule in solution is random, so is the orientation of the moments, and consequently, so is the orientation of the tensorsâ,b.
The macroscopic polarization P of the chiral isotropic liquid takes the form P ¼ Nhmi; (20) where N is the number of molecules and the brackets h.i denote the isotropic average (uniform over all directions). Thus, in order to calculate P, we need to consider the isotropic average It can be shown 51,52 that the isotropic average of any secondorder tensorû equals Now, to obtain the characteristic macroscopic propertymacroscopic susceptibility tensorĉ, we can make use of the standard denition This equation can be solved for monochromatic plane waves E = E 0 e −ik$r , where E 0 is the complex electric eld amplitude, k is the wavevector and r is the position vector. Then, rot E h V × E = −ik × E. The cross product can be elegantly avoided by rewriting the vector k as follows, ik Â EhiK$E; whereK ¼ 2 6 6 4 0 Àk z k y k z 0 Àk x Àk y k x 0 is the wave vector tensor.
The macroscopic susceptibilityĉ, which is a macroscopic experimental observable tied to the refractive index of the medium vian 2 =3 = 3 0 (1 +ĉ), then holds respective elements will not directly correspond to the above quantities but will be their measures, each separately lying on the matrix secondary diagonal. Comparing the matrix (29) with (15)- (17), the elements of interest correspond to those that work exclusively with circular polarizations. We offer an illustrative example regarding chiral solution measurements: If we use an ideal transparent cuvette without residual stress in the faces, the response of the raw experimental Mueller matrix will arise only on the secondary diagonal. Elsewise, the material stress in the cuvette (induced linear birefringence in terms of polarimetry) will give rise to signals in the elements m 24 , m 34 , m 42 , and m 43 . We illustrate this most complex case of matrix (29) in Fig. 4b. The drawn example would correspond, for example, to the response of a chiral compound with the exciton coupling causing negative bisignate CD 17 and Kramers-Kronig related CB.
Even if we do not have a physical model, from a general acquaintance with the behavior of the polarimetric response in the Mueller matrix, it is already possible to determine the elements of our interest. Since the individual polarization effects do not mix up in most Mueller matrices, M diff may stand for a solid clue for the experimentalist.
Note that the presented approach is rather illustrative and should be applied only to the discussed specic cases, importantly, certainly not to the reection Mueller matrices.

Conflicts of interest
There are no conicts to declare.